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THE  USE  OF  MEASUREMENTS  IN  PREDICTING  THE 


PERFORMANCE  OF  TROPOSPHERIC  COMMUNICATION  CIRCUITS 


by 


A.  P.  Barsis,  P.  L.  Rice  and  K.  A.  Norton 


Summary 

The  uncertainty  in  predicting  the  performance  of  tropospheric 
communication  circuits  may  be  expressed  by  a  variance,  assuming  a 
log-normal  distribution  of  the  differences  between  calculated  and 
measured  values  of  field  strength  or  basic  transmission  loss,  if  a 
number  of  paths  with  identical  geometrical  and  other  prediction  param¬ 
eters  are  considered. 

It  is  shown  to  what  extent  this  variance  may  be  reduced  by  taking 
into  account  transmission  loss  measurements  over  the  same  path  for 
which  the  prediction  is  intended,  or  over  a  similar  path.  The  "improved 
prediction"  is  obtained  by  assigning  weights  to  calculated  and  observed 
values,  respectively,  and  forming  a  weighted  average.  The  method  is 
demonstrated  by  examples,  and  the  results  are  used  to  calculate  new 
values  of  service  probability  on  the  basis  of  the  improved  prediction. 
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Comparison  of  service  probabilities  for  predictions  based  on  the  cal¬ 
culation  above,  and  for  predictions  taking  into  account  the  results  of 
transmission  loss  measurements,  enables  the  designer  to  judge 
whether  the  cost  of  a  measurement  program  is  justified  in  terms  of 
possible  savings  in  equipment  for  a  final  installation. 

1.  Introduction 

In  a  recent  report  by  the  authors-^/a  method  was  given  to 
predict  the  performance  of  a  tropospheric  communication  system  by 
the  use  of  semi -empirical  methods,  and  to  evaluate  the  probability  of 
its  success  by  the  introduction  of  the  concept  of  service  probability. 
Briefly,  service  probability  has  been  defined  as  the  probability  of  ob¬ 
taining  a  specified  grade  of  service  or  better  during  a  given  percentage 
of  the  hours  in  a  specified  period  of  time,  usually  taken  to  be  all  hours 
of  the  year.  The  service  probability  is  designated  as  F(t),  and  the 
percentage  of  the  hours  during  the  period  considered  is  termed  time 
availability  and  designated  by  p.  The  specified  grade  of  service  is 
predetermined  in  terms  of  an  hourly  median  r.m.s.  carrier  to  r.m.  s. 
noise  ratio  which  takes  into  account  short-term  signal  fading  and  the 
modulation  capabilities  of  the  circuit.  All  three  concepts,  namely, 
grade  of  service,  time  availability,  and  service  probability  are  brought 

together  in  a  single  graphical  representation.  The  final  graph  has  two 
\J  A.  P.  Barsis,  K.  A.  Norton,  and  P.  L.  Rice,  "Predicting  the  per- 

formance  of  long  distance  communication  circuits,  "  NBS  Report  60  32, 

December  1958 
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probability  scales  with  the  ordinate  showing  the  time  availability  p  in 
percent  of  all  hours  and-  the  abscissa  showing  the  service  probability 
F(t).  Each  curve  on  the  graph  corresponds  to  a  system  with  fixed 
operating  characteristics  such  as  transmitter  power,  antenna  size,  etc.  , 
and  may  therefore  be  labelled  by  a  figure  representing  installation  and 
operating  costs.  Various  systems  may  thus  be  compared  on  the  basis 
of  cost  versus  service  probability.  In  other  words,  the  cost  of  various 
systems  may  be  compared  using  the  probability  of  success  for  each 
system  as  a  yardstick. 

From  a  study  of  the  methods  given  in  Ref.  1  and  from  an 
inspection  of  service  probability  displays,  it  is  shown  that  there  is  an 
allowance  made  for  uncertainties  in  the  prediction  process.  This  un¬ 
certainty  is  expressed  in  terms  of  a  normal  distribution  with  standard 
deviation  °'rc(P)‘  The  standard  deviation,  crrc(p),  includes  a  term 

a  (p)  for  errors  in  the  estimate  of  transmission  loss,  as  well  as  a 
c 

term  crr  which  allows  for  variations  in  equipment  performance.  Fig.  1 
shows  these  uncertainties  as  derived  in  Ref.  1.  If  cr  (p)  can  be  reduced 
below  the  values  given,  then  the  estimate  of  the  power  required  for  high 
service  probabilities,  such  as  F(t)  =  0.95,  will  be  reduced  accordingly. 
Also,  service  probability  curves  determined  by  these  methods  will  be 
less  steep,  and  the  expected  value  p  for  F(t)  =0.5  would  constitute  a 
much  better  estimate  since  it  would  not  differ  as  much  from  the  esti- 
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mates  for  higher  and  lower  values  of  service  probability.  The  question 

is  how  much  the  standard  deviation  crrc(p)  can  be  decreased  by  taking 

into  account  the  results  of  transmission  loss  measurements  made  for 

various  periods  of  time.  The  uncertainty  in  the  assumption  of  median 

r.m.s.  carrier  to  median  r.m.  s.  noise  ratio  R  and  of  the  effective  noise 

figure  F,  which  has  been  expressed  by  the  standard  deviation  <r ,  will 

not  be  affected  by  the  availability  of  transmission  loss  or  field  strength 

measurements,  and  only  the  crc(p)  portion  of  the  standard  deviation 

o’  (p)  can  be  reduced.  <r  (p)  as  a  function  of  p  is  shown  by  curves  (a) 
i  c  c 

and  (c)  on  Fig.  1,  together  with  the  corresponding  weights.  The  use  of 
the  weights  will  be  discussed  later. 

2.  Statistical  Basis  for  Prediction  Improvement 

Analogous  to  the  prediction  uncertainty  v^(p)  which  was  assigned 
to  calculated  values  of  field  strength,  a  standard  deviation  <rQe(p)  is 
assigned  to  observed  field  strength  values.  An  improved  prediction  is 
then  obtained  by  computing  weighted  averages  of  predicted  and  measured 
values.  In  conformance  with  statistical  practice,  weights  are  defined 
as  the  reciprocals  of  the  variances.  The  variance  cr£(p)  has  already 
been  defined;  cr  (p)  and  the  reciprocal  of  cr^(p)  which  is  the  weight  w  (p) 
are  shown  on  Fig.  1.  An  expression  for  <r0g(p)  will  be  derived  in  the 
next  section.  Most  of  the  terms  and  symbols  used  are  tabulated  in 


Appendix  I. 
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The  general  problem  considered  here  is  the  improvement  of 
predictions  made  for  a  specific  tropospheric  propagation  path,  using 
measurements  obtained  over  the  same  or  a  slightly  different  path.  It 
will  often  be  the  case  that  measurements  are  made  only  during  several 
weeks  or  months,  while  the  improved  prediction  should  be  applicable  to 
all  hours  of  the  year. 

As  the  calculations  for  tropospheric  communication  systems 
are  usually  in  terms  of  basic  transmission  loss  the  method  of 

prediction  improvement  will  also  be  treated  in  these  terms.  It  should 
be  remembered,  however,  that  the  two  quantities,  field  strength  and 
basic  transmission  loss,  are  directly  related,  and  an  hourly  median 
value  of  field  strength  E(p)  exceeded  during  at  least  p %  of  all  hours 
is  equivalent  to  a  basic  transmission  loss  value  L^P)  not  exceeded  during 
at  least  p %  of  all  hours.  In  the  following  discussion,  all  L  values  are 
hourly  medians,  and  the  p-dependence  is  implied,  although  the  symbol  (p) 
will  be  dropped  for  convenience. 

If  an  observed  transmission  loss  value  L  and  a  calculated 


oa 

value  L  ,  both  in  decibels,  are  available  for  Path  a,  the  improved 

C  cL 

prediction  Lta  is  given  by: 


Jta 


woa  Loa  +  wca  Lca 

w  +  w 
oa  ca 


(1) 


where  w„„  and  w 

oa  ca 


are  weights  assigned  to  the  observed  and  calculated 
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values  for  Path  a. 

Suppose,  however,  that  an  improved  prediction  is  desired 

for  Path  b,  for  which  measured  values  are  not  available.  Calculated 
values  Lca  and  L  ^  are  available  for  both  paths,  and  we  may  estimate 
as  follows: 


Ltb  =  Lta+<Lcb  -Lca>  <2> 

The  measured  value  L,  „  is  the  sum  of  a  "true  value"  L 

oa  a 

(independent  of  p)  and  a  component  5  of  the  measurement  error 

cL 

associated  with  the  fact  that  the  measurements  are  made  over  a  finite 
period  of  time* 


L 


oa 


L  +6. 


(3) 


The  computed  values,  L  and  L  ,  are  the  sum  of  true  values 

ca  cb 

L  and  L,  and  prediction  errors  e  „  and  eu: 
a  b  a  b 


L  =  L  +  e  , 
ca  a  a 


L 


cb 


=  L,  +  € 


b 


(4) 


The  expected  values  of  6a>  and  are  zero,  the  expected 

2  2 
variance  of  L  is  cr  =  l/w  the  expected  variance  of  L  K  is  tr  u, 
ca  ca  7  ca  r  cb  cb 

*  In  addition  to  this  measurement  error  there  will  always  be  some  un¬ 
certainty  in  the  accuracy  of  the  measurements  due  to  calibration 
errors;  this  latter  uncertainty  is  taken  to  have  a  variance  cr^  and  will 
be  introduced  later. 
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and  the  expected  variance  of  L  is  <r  ^  =  l/w 

r  oa  oa  •  oa 

It  will  be  assumed  that  the  correlation  between  L  and  either 

oa 

L  or  L  .  is  zero.  The  correlation  between  L  and  L  ,  will  be 
ca  cb  ca  cb 

denoted  by  pa^.  As  a  practical  matter,  it  will  sometimes  be  convenient 

to  set  pak  equal  to  either  zero  or  one,  depending  on  whether  the  paths 

are  very  different  or  almost  alike. 

With  the  above  assumptions,  the  variance  about  the  estimate  of 

Lj-b  in  (2)  may  be  obtained.  It  is  shown  in  Appendix  II  that  the  corre- 

lation  between  L.„  and  L  is  cr  /<r  ,  where  oy  =  l/(w  +  w  )  is  the 

ta  ca  ta'  ca  ta  M  oa  ca' 

variance  of  the  estimate  .  This  is  the  estimate  used  when  measure - 

ta 

ments  are  made  over  the  path  for  which  predictions  are  desired,  and 
the  remaining  sections  of  this  paper  will  be  devoted  to  a  derivation  of 
crt^  .  The  correlation  between  Lta  and  Lc^  is  (0‘ta/crca)pab“  Then  the 
variance  associated  with  the  improved  prediction  (2)  for  Path  b  is 
given  by 

“tb  =  "cb  +  °c i  -  °'ta  -  2'rca°-cbPab(1  ‘  <rtJ/<rCi)  (5) 


The  detailed  derivation  of  (5)  is  given  in  Appendix  II. 

If  Path  a  is  the  same  as  Path  b,  then  =  Lta„  the  variances 

=  1  !  in  thi  s  raR 

ab 


<r_| f  and  o’  ^  are  equal,  and  p  =  1;  in  this  case,  (5)  reduces  to  cr^ 


ca 


<r  as  it  should, 
ta 


It  is  shown  in  Appendix  II  that  the  correlation  coefficient  pa^ 
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between  estimates  L  _  and  L  ,  over  the  two  Paths  a  and  b  has  to  be 

ca  cb 

at  least  0.  5  in  order  to  reduce  the  standard  error  of  estimating  Ltb 
below  the  standard  error  of  L  ^  In  general,  it  will  be  quite  difficult 
to  make  an  accurate  estimate  of  the  correlation  coefficient  p  .  .  The 
results  of  the  study  in  Appendix  II  show  that  pa^  has  to  be  very  close 
to  unity  in  order  to  obtain  a  substantial  improvement  in  the  estimate  of 
L'tb  over  the  estimate  obtained  without  any  measurements.  Avail¬ 

able  data  which  may  be  used  to  estimate  p  are  contained  in  papers  by 
Kirby  and  Capps—/  and  Kirby.1/.  It  is  shown  there  that  for  a  fixed  trans¬ 
mitter  location  the  serial  correlation  coefficient  between  field  strength 

values  observed  at  adjacent  measuring  locations  drops  to  0.  5  within  a 

4/ 

distance  of  one -half  to  one  mile.  On  the  other  hand,  Kirby,  et  al—', 
have  shown  that  there  is  relatively  high  correlation  between  different 
frequencies  received  over  essentially  the  same  paths.  All  these  results 


were  obtained  on  much  shorter  paths  than  are  normally  used  for 

tropospheric  scatter  communication  systems,  and  may  not  be  directly 

applicable.  To  the  extent  that  these  measurements  are  applicable  it 

would  appear  that  path  loss  measurements  made  on  paths  with  one  or 

the  other  of  their  terminals  separated  by  as  much  as  one  mile  from  the 
2/  R.  S.  Kirby  and  F.  M.  Capps,  "Correlation  in  VHF  propagation  over 

irregular  terrain,"  IRE  Trans.  PGAP,  Vol.  AP-4,  No.  1,  pp.  77-85, 
January  1956 

3/  R.  S.  Kirby,  "Measurement  of  service  area,  for  television  broad¬ 
casting,  "  IRE  Trans.  PGBTS-7,  pp.  23-30,  February  1957 
4 /  R.  S  .  Kirby,  H.  T.  Dougherty,  and  P.  L.  McQuate,  "VHF  propaga¬ 
tion  measurements  in  the  Rocky  Mountain  region,  "  ERE  Trans. 
PGVC-6,  pp.  13-19,  July  1956  (see  Fig.  6) 
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terminals  of  the  path  for  which  a  prediction  is  required  will  be  of  little 
or  no  direct  value. 

The  method  of  obtaining  the  variance  ^^(P)  and.  weight  w^(p) 

for  a  set  of  measurements  and  the  procedure  for  applying  them  to  all 

day-all  year  estimates  will  be  discussed  in  the  following  section. 

3.  Estimating  the  Standard  Deviation  croe(p) 

The  variance  0"o|(p)  associated  with  an  observed  field  strength 

value  EQ(p),  or  the  equivalent  basic  transmission  loss  L^Q(p)  may  be 

2  2 

separated  conveniently  into  two  parts  <r  (p),  and  cr  ,  so  that 

O  c 


O'  ?(p)  =  o-* 2(p)  +  fl¬ 


oe 


(6) 


In  this  representation  cre  is  the  "equipment,  "  or  reading  error, 
taken  to  be  independent  of  the  percentage  of  hours,  p.  The  part  cr0(p), 
however,  is  considered  to  be  a  function  of  the  time  availability  p  in 
percent  of  hours,  the  number  of  days  of  measurement,  n,  and  the 
angular  path  distance,  0.  Thus  cr^(p)  is  that  part  of  the  variance  of 
observed  values  Lj-^P)  independent  of  the  error  tre  introduced  by  the 
measuring  equipment  and  procedures. 

The  weight  wq  =  l/cro(p)  associated  with  measurements  extending 
over  a  short  period  of  time  is  less  than  the  weight  associated  with  the 
median  of  basic  transmission  loss  values  observed  over  a  long  period. 

If  only  one  hourly  median  were  available  for  some  path,  for  instance, 
it  would  be  the  only  available  estimate  of  L^o(50).  The  variance 
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ctq(P  =  50%,  N  =  1  hour)  associated  with  one  hourly  median  is  simply 
the  variance  of  all  hourly  medians  that  might  have  been  recorded  over 
the  path. 

Analysis  of  a  large  amount  of  tropospheric  propagation  data  has 
resulted  in  the  determination  of  the  dependence  of  cro(50)  on  the  number 
of  days  of  observation.  There  is  so  much  more  correlation  between 
measurements  made  during  successive  hours  of  one  day  than  there  is 
between  measurements  made  at  the  same  hour  on  successive  days, 
that  the  total  number  of  days  is  a  better  measure  of  sample  size  than 
the  total  number  of  hours.  One  might  conclude  then  that  one  hour  of 
observation  a  day  for  a  great  many  days  is  worth  almost  as  much  as  a 
continuous  period  of  recording.  However,  an  exception  to  the  above 
statement  needs  to  be  made  in  order  to  take  care  of  the  commonly 
observed  diurnal  variation  of  radio  transmission  loss.  As  in  earlier 
work.— 1-ilA/,  the  hours  of  the  year  are  grouped  into  "time  blocks" 
which  divide  the  day  into  four  parts  and  the  year  into  two  parts,  winter 

and  summer.  An  examination  of  the  dependence  of  0^(50)  on  the  number 

5/  P.  L.  Rice,  A.  G.  Longley,  and  K.  A.  Norton,  "Prediction  of  the 
cumulative  distribution  with  time  of  ground  wave  and  tropospheric 
wave  transmission  loss,  "  NBS  Report  5582,  June  1958 
6 j  K.  A.  Norton,  P.  L.  Rice,  L.  E.  Vogler,  "The  use  of  angular 
distance  in  estimating  transmission  loss  and  fading  range  for 
propagation  through  a  turbulent  atmosphere  over  irregular  terrain,  " 
Proc.  IRE,  43,  No.  10,  pp.  1488  -  1526,  October  1955 
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of  days,  n,  was  then  made  for  each  time  block.  Table  I  below  lists 
all  time  blocks. 


Table  I 

Time  Blocks 

No. 

Month 

Hours 

1 

November  -April 

0600-1300 

2 

November  -April 

1300-1800 

3 

November  -April 

1800-2400 

4 

May-October 

0600-1300 

5 

May-October 

1300-1800 

6 

May-October 

1800-2400 

7 

May-Oc  tober 

0000-0600 

8 

November  -April 

0000-0600 

It  has  also  been  found  convenient  to  combine  the  four  time 
blocks  pertaining  to  the  winter  season  (1,  2,  3  and  8)  and  designate 
this  combination  as  a  time  block  group.  Similarly,  the  combined 
Time  Blocks  4,  5,  6  and  7  form  another  time  block  group  pertaining 
to  the  summer  season.  From  the  analysis  of  a  large  amount  of  propa¬ 
gation  data  more  fully  described  in  Ref.  6,  an  empirical  relation 
has  been  determined  between  the  hourly  median  value  of  field  strength 
(or  basic  transmission  loss)  calculated  for  Time  Block  2  (the  average 
winter  afternoon  hour)  and  the  hourly  median  applicable  to  any  de¬ 
sired  percentage  of  hours  within  a  specific  time  block  or  time  block 
group. 

This  relation  is  given  by  the  V(p,  0)  curves  shown  on  Figs.  2, 
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3  and  4,  where  the  hourly  median  basic  transmission  loss  L^m(p) 
corresponding  to  the  field  strength  value  exceeded  during  p%  of  all 
hours  is  given  by 

Lbm'P>  =Lbm<50>  -  V<P’9>  <7> 

In  (7),  L^m(50)  is  always  the  Time  Block  2  median  value  (calculated 
in  accordance  with  the  methods  given  in  Ref.  5),  and  the  value  V(p,  9) 
is  read  from  Figs.  2,  3  or  4,  depending  on  whether  the  distribution 
of  hourly  medians  is  desired  for  the  winter  time  block  group  (later 
designated  by  "F"),  the  summer  time  block  group  (similarly  designated 
by  "S").  or  for  all  hours  of  the  year  (designated  by  "A"). 

Graphs  similar  to  Figs.  2,  3  and  4,  but  for  individual  time 
blocks  are  contained  in  Ref.  5. 

From  an  analysis  of  long-term  data  from  12  propagation  paths 
having  a  mean  angular  distance  of  9  =  18.  2  milliradians ,  a  variance 
curve  of  cr00(59)  versus  n  for  9  =  18.2  milliradians  was  obtained  and 
plotted  on  Fig.  5,.  That  is,  each  long  recording  period  was  broken  up 
into  several  shorter  periods  of  various  lengths,  and  the  variance 
among  medians  of  groups  of  hourly  medians  was  determined  for  each 
of  several  sample  sizes.  As  an  example,  for  n  =  36  days,  Fig.  5 
reads  croo(59)  =  19.  6  db  ;  this  is  the  variance  of  the  medians  of  five 


groups  of  189  hourly  medians  each,  representing  five  36 -day  recording 
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periods.  Actually,  the  figure  shows  averages  of  values  obtained  for 
a  number  of  paths  in  Time  Block  2,  winter  afternoons.  The  shape  of 
the  variance  curve  for  other  time  blocks  was  not  found  to  be  signifi- 
cantly  different.  The  level  of  croo(50)  for  n  =  0.  2  was  determined  from 
the  estimates  of  long-term  variability  for  all  hours  as  given  in  Fig.  4, 
assuming  that  n  =  0.  2  corresponds  to  one  hour,  as  it  does  for  Time 
Block  2.  Fig.  5  also  contains  a  curve  for  the  weight  woo(50),  the 
reciprocal  of  cro^(50),  as  a  function  of  n.  This  weight  curve  is  added 
for  convenience,  if  needed  directly  in  calculations. 

Available  data  indicate  that  the  variances  o-£(50)  are  much  more 
a  function  of  angular  distance  than  the  prediction  uncertainty 
discussed  in  Sec.  6  of  Ref.  1,  and  shown  here  on  Fig.  1.  Curves  of 
<Tq(50)  versus  n  for  various  angular  distances  are  roughly  parallel, 
and  their  mutual  relation  corresponds  approximately  to  the  V(p,  9) 
curves  shown  on  Figs.  2,  3  and  4.  A  factor  f(9)  has  been  determined 
on  the  basis  of  the  all  day-all  year  V(p,  9)  curves,  by  which  the 
variance  (roo(50)  shown  on  Fig.  5  corresponding  to  an  angular  distance 
of  18.  2  milliradians  has  to  be  multiplied  to  obtain  the  variance 
cr^(59)  for  a  given  angular  distance  9.  The  factor  f(9)  as  a  function  of 
9  is  shown  on  Fig.  6,  and  the  relation  may  be  expressed  by  the 


equation 
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%<50>  =<roo<50>f<9>  (8) 

The  following  paragraphs  explain  how  the  variance  tr^p)  is 
obtained  for  values  of  p  other  than  50%.  Here,  cr^(p)  is  the  variance 
of  observed  hourly  median  field  strength  values  EQ(p)  exceeded  for 
p%  of  all  hours.*  On  the  average,  EQ(p)  and  Lbo(p)  are  related  to 
Eq(  50)  and  Ebo(50)  as  shown  in  (  7)  above: 

Eq(p)  =  Eq(50)  +  AV  (p,  0) ,  (9) 

or 

Lbo(P)  ^  Lbo(50)  "  AV(p,0)  (9a) 

where 

AV(p,  9)  =  V(p,  0)  -  V( 50 ,  0)  (10) 

It  should  be  noted  that  here  Eq( 50)  or  Lbo  (50)  refers  to  the 
median  of  the  time  block  or  time  block  group  considered.  It  is  also 
assumed  that  errors  in  Eq(  50)  and  AV(p,  0)  are  uncorrelated,  and  that 
the  variance  <^y  (P)  is  proportional  to  (AV)^: 

<rAv(P>  =k2(AV)2  (11) 

It  follows  that  cr^(p)  may  be  written 

*Again,  a  field  strength  value  E(p)  exceeded  during  p%  of  all  hours  is 

equivalent  to  a  transmission  loss  value  L(p)  not  exceeded  during  the 
same  p%  of  all  hours. 
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<Tq(P)  =  <^(50)  +  k2(AV)2  (12) 

where  cr^(50)  is  obtained  from  (8),  using  Figs.  5  and  6,  while  AV  is 
obtained  from  Fig.  2,  3  or  4,  whichever  is  applicable  to  the  measure¬ 
ment  period  considered. 

The  reasonableness  of  the  above  assumptions  was  verified  and  a 
value  of  0.09  for  was  determined  using  the  same  Ohio  data  dis¬ 
cussed  in  the  appendix  to  Ref.  1.  Certain  parameters  for  these  data 
are  listed  in  Table  II,  which  gives  values  of  (Tq(50)  determined  from 
(8)  for  each  of  two  frequencies  and  two  distance  ranges. 

The  number  of  days  of  measurement  shown  in  Col.  6  of  Table  II 
represents  a  full  year  of  random  sampling  at  12  5  miles  and  seven 
months  of  random  sampling  at  8  5  miles. 

Observed  values  of  o-^(p)  for  each  set  of  data  were  normalized 

o 

to  agree  at  p  =  50%  with  the  theoretical  value  listed  in  Table  II  by 
adding  or  subtracting  an  appropriate  variance.  These  normalized 
variances  o'Q(p)  for  the  four  sets  of  data  were  then  averaged;  these 
averages  are  listed  in  Table  III.  The  V(p,  9)  values  in  Table  III  are 
based  on  the  average  9  =  2  5.  7  milliradians  for  all  paths. 


Data  for  Determination  of  cTq(50) 
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Table  III 


P 


1% 

10% 

50% 

90% 

99% 

ro(p> 

37.  59 

19.  61 

12.  74 

21.  22 

29.  54 

V(p,  0) 

22.  2 

13.  0 

3.  8 

-  4.0 

-10.  2 

(from  Fig. 

4) 

AV 

(from  (10)) 

18.  4 

9.  2 

0 

-  7.8 

-14.  0 

k2 

(from  (12)) 

0. 0734 

0. 0812 

- 

0. 1394 

0.0856 

k 

0.  271 

0.  285 

- 

0.  373 

0.  293 

p(P) 

(from  ( 1  3) 

-  0.  18 

-  0.06 

- 

+  0.  16 

-  0.06 

below) 

The  average  value  of  in  Table  III  is  approximately  0.09,  as 
determined  by  assuming  no  correlation  between  errors  in  Eq(50)  and 
AV.  If  the  value  of  k  is  fixed  at  0.  3  and  it  is  assumed  that  cr  .  -r  =  k  I  AV 
then  the  calculated  correlation  coefficient  p(p),  using  the  data  of 
Table  III,  is  seen  to  be  essentially  zero  at  each  percentage  level.  For 
this  calculation,  the  relation 

„2(p)  =  0-2(50)  +  <rA^(p)  +  2p(p)<ro(50)<rAV(p)  (13) 

was  used. 

2 

Thus  (12)  is  justified,  with  0.09  as  the  value  for  k  . 
Recapitulating,  we  now  have  an  estimate  of  the  variance  cr2(p)  to 
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be  assigned  to  a  measured  field  strength  value  E^(p)  exceeded  by  p% 
of  all  hourly  medians.  Combining  (9),  (10)  and  (12): 

o-o(P)  =  ^oo<50>£(9>  +  °-°9[V(p,  6)  -  V(50,  0)]  (14) 

where  Figs.  5  and  6  and  one  of  the  appropriate  curves  of  Figs.  2,  3 

2 

or  4  are  used  to  evaluate  (Tq(p)  from  (14). 

The  preceding  paragraphs  in  this  section  assume  that  measure¬ 
ments  are  made  in  the  same  time  block  for  which  estimates  are  de¬ 
sired.  It  will  usually  be  the  case  that  measurements  will  be  made  for 
a  period  of  weeks  or  months  either  in  summer  or  winter,  and  that 
estimates  will  be  desired  for  all  year.  A  generalization  of  (12)  is 
here  made  which  depends  upon  a  generalization  of  (9)  or  (9a)  inherent 
in  prediction  methods  described  in  Ref.  5.  Suppose  we  have  available 
a  measured  cumulative  distribution  for  the  summertime  and  desire 
estimates  of  L(p)  for  all  year.  Denote  by  Vg(p,  9)  the  function  V(p,  0) 
plotted  in  Fig.  3  for  the  summer  months  May  through  October,  and 
denote  by  V^(p,  9)  the  function  V(p,  9)  plotted  for  all  hours  in  Fig.  4. 
Then  re -define  AV  as  follows;  with  the  double  subscript  denoting  that 
it  is  the  difference  between  a  V(p,  9)  value,  V^( p,  0)  in  time  block 
group  A  (all  year),  and  a  V(p,  9)  value,  Vg(p,  0)  in  time  block  group  S 
(summer): 


AVAS  2  AVAS(P’9)  =  VA(P’0)  -  vs(p,  6) 


(15) 
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An  estimate  of  L^^(p)  for  all  hours  is  obtained  from  the  mea¬ 


sured  value  of  Lj-^P)  using  the  following  relation: 


LbA,P>  =  LbS(P)  -  AVAS(P'  S) 


(16) 


This  is  the  analog  to  (8b),  and  (10)  may  be  rewritten  for  the  summer 
period: 


(10a) 


AVS  =  AVs(p,  9)  -  Vs(p,  9)  -  Vs(50,  0) 


Our  extension  of  (12)  is  then 


°o(P)  55  °AS(P)  =  "o(5°)  +  kZ  [  (AVS)2  +  (AYAS)2] 

=  cr Qq( 50)f ( 0)  +  0.  09  {[  Vs(p,  0)  -  Vg(  50 ,  0)]  2 
+  [VA(p,  0)  -  Vg(p,  0)]  2) 


(17) 


Equation  (17)  now  replaces  (14),  which  it  generalizes,  so  that 


a  measurement  L.  (p)  made  in  the  summertime.  Equation  (16)  is 
bS 

used  to  estimate  a  value  for  all  year  from  measurements  made  during 
part  of  the  year.  If  measurements  are  made  in  the  winter,  change  the 
subscript  S  to  F  in  the  preceding  equations  and  use  Fig.  2  instead  of  3. 
In  cases  where  measurements  are  made  during  individual  time  blocks 
only,  the  appropriate  V(p,  0)  curves  from  Ref.  5  have  to  be  used  in¬ 


stead  of  Figs.  2,  3  or  4. 
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4.  Application  of  the  Standard  Deviation  croe(p)  to  Predictions ,  and 
Prediction  Improvement 

The  standard  deviations  aQ(p)  calculated  in  accordance  with  the 
methods  given  in  the  preceding  section  characterize  the  prediction 
uncertainty  associated  with  any  observed  field  strength  or  trans¬ 
mission  loss  value  if  the  observed  value  is  not  subject  to  equipment  or 
reading  errors.  Section  3  introduced  the  symbol  <xe  to  characterize 
such  errors,  considered  to  be  independent  of  the  measured  field 
strength  or  transmission  loss  value.  Under  the  reasonable  assumption 
that  the  prediction  uncertainty  and  the  equipment  errors  are  un¬ 
correlated,  the  weight  wQ(p)  assigned  to  an  observed  value  L^^p)  is 
is  given  by 


1 


(6a) 


WQ(P) 


For  the  usual  transmission  loss  measurement  and  analysis 
procedures,  a  value  of  2  decibels  constitutes  a  good  estimate  of  c^. 

As  the  standard  deviation  &oe(p)  now  characterizes  the  un¬ 
certainty  associated  with  the  observed  value  L  (p)  for  each  per- 

bo 

centage  point  p,  it  may  be  used  directly  to  draw  a  service  probability 
graph  in  accordance  with  the  method  described  in  Ref.  1«  This  means 
that  the  performance  of  a  system  may  be  estimated  on  the  basis  of 


measurements  alone  without  any  calculation  of  basic  transmission  loss 
at  all.  For  this  procedure  it  is  only  necessary  to  modify  croe(p)  by 
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taking  into  account  the  uncertainty  in  estimating  system  equipment 
parameters  just  as  was  done  in  the  theoretical  prediction  case  (see 
p.  10  of  Ref.  1).  The  total  uncertainty  cror(p)  in  this  case  is  given  by 
an  equation  analogous  to  (4)  of  Ref.  1: 


(18) 


where  <r  ,  usually  taken  as  2  decibels,  characterizes  the  uncertainty 
in  estimating  system  equipment  parameters.  Service  probabilities 
may  now  be  calculated  in  accordance  with  (5)  and  (6)  of  Ref.  1. 
Observed  transmission  loss  values  L^Q(p)  are  related  to  system  power 
requirements  by  appropriate  constants  depending  on  the  grade  of 
service  required  if  all  quantities  are  expressed  in  decibels.  A  con¬ 
venient  expression  for  this  is  given  by  Norton-^/, 


Pto  =  Lbo  -  Gp  +  F  +  Lt  +  R  +  B  -  204 


(19) 


where 


=  transmitter  power  in  decibels  relative  to 
one  watt 


L^q  =  basic  transmission  loss  over  the  path  in  decibels 


G  =  path  antenna  gain  in  decibels 


F  =  effective  receiver  noise  figure  in  decibels 


7 /  K.  A.  Norton,  "Point-to-point  radio  relaying  via  the  scatter  mode 

of  tropospheric  propagation,  "  Trans.  IRE  on  Comm.  Systems, 
Vol.  CS-4,  No.  1,  pp.  39-49,  March  1956 
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Lj.  =  transmission  line  and  associated  losses,  in  decibels 

R  =  pre -detection  r.  m.  s„  carrier  to  r.m.s.  noise  ratio  in 

decibels,  for  the  bandwidth,  type  of  modulation,  and 
order  of  diversity  used 

B  =10  log 1Qb,  where  b  is  the  total  bandwidth  in  cycles  per 
second  (including  the  effect  of  drift  between  transmitter 
and  receiver  oscillators) 

204  =  a  constant  equal  to  -10  log^kT,  where  k  is  Boltzmann's 
constant,  and  the  reference  temperature  T  is  taken  to  be 
288.  44°  Kelvin 

The  usual  procedure  is  to  utilize  measurements  in  conjunction 
with  predictions,  in  order  to  reduce  the  prediction  uncertainty.  Re¬ 
capitulating  the  results  obtained  in  Section  2  above,  it  is  necessary  to 


distinguish  cases  where  measurements  are  available  for  exactly  the 
same  path  for  which  the  prediction  is  desired  from  the  ones  where  the 
paths  are  similar,  or  in  the  same  general  area. 

In  the  first  case,  the  improved  prediction  is  given  by  (1).  It  may 
be  re-written  in  terms  of  basic  transmission  loss  values: 


w0(p)Lb0  +  w  (p)Lbc 

L  (p)  = - { ta.) 

bt  wQ(p)  +  wc(p) 

The  weight  and  the  variance  assigned  to  the  improved  prediction 
Lbt<P)  is  §iven  bY 


wt(p)  -  — -  =  wQ(p)  +  wp(p) 

Tt<P) 


(20) 


with  w  (p)  shown  on  Fig.  1  as  a  function  of  p. 
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o"t(p)  is  a  measure  of  the  remaining  prediction  uncertainty  characteris¬ 
tic  of  the  improved  basic  transmission  loss  values  L^p)  just  as 
cr  (p)  was  a  measure  of  the  original  prediction  uncertainty  pertaining 
to  the  theoretical  values  L^c(p).  If  system  power  requirements  are 
considered,  the  term  crr  has  to  be  added  to  cr^.  (p)  in  order  to  obtain 
a  measure  of  the  remaining  total  prediction  uncertainty  taking  into 
account  the  uncertainty  in  estimating  the  system  parameters  F  and  R 
as  discussed  on  p.  10  of  Ref.  1.  Thus,  the  remaining  prediction  un¬ 
certainty  is  characterized  by  a  standard  deviation  crt(p)  as  a  function 
of  p,  if  hourly  medians  of  transmission  loss  or  field  strength  values 
are  considered,  or  by  a  standard  deviation  crtr(p),  if  predictions  are 
to  be  made  for  the  performance  of  a  complete  system.  o"^r(p)  is 
given  by 

o^(p)  =  o-f(P)  +  (21) 

For  the  second  case,  where  measurements  are  available  for 

one  path  (Path  a),  whereas  a  prediction  is  desired  for  another  path 

(Path  b),  the  improved  prediction  on  the  basis  of  measurement  for 

Path  a  and  its  variance  is  first  determined  in  accordance  with  (la)  and 

(20)  above.  To  this  improved  prediction  L^ta»  the  difference  of  the 

calculated  values  L.  and  L,  ,  for  the  two  paths  is  added  as  in  (2) 

bca  bcb 


in  Sec.  2  above; 
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L,  ,  -  L,  +  (L  -  L,  i  ) 
btb  bta  bca  bcb' 


(2a) 


The  variance  cr^  characterizing  Lbj_b  is  now  obtained  in 
accordance  with  (5)  using  an  appropriate  assumption  for  the  correla¬ 
tion  coefficient  pa^«  For  the  design  of  a  system,  the  variance  c has 
to  be  added  to  a tb,  as  in  (21). 


crtrb(P)  =  ^tb^P)  +  *r 


(21a) 


Having  obtained  Lbt(p)  or  Lbt±)(p)  and  crtr(p)  or  <rtrb(p)  as 
functions  of  p,  service  probability  values  may  be  calculated  in  the 
same  way  as  was  done  for  the  theoretical  prediction  described  in 
Ref.  1. 

The  procedure  will  be  demonstrated  by  examples  in  the  next 
section. 

5.  Examples 

In  this  section  the  prediction  improvement  method  is  illustrated 
by  the  use  of  two  examples.  In  the  first  example,  measurements 
made  over  a  200 -mile  path  between  North  Africa  and  Spain  are  used  to 
demonstrate  prediction  improvement.  This  is  the  sample  path  dis¬ 
cussed  in  Ref.  1.  The  second  example  illustrates  the  effect  of  in¬ 
creasing  the  length  of  the  measurement  period,  and  a  portion  of  long¬ 
term  measurement  results  (data  obtained  from  an  FM  broadcasting 
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station  in  the  U.  S.  )  is  used. 

(a)  Prediction  Improvement  for  the  Sidi  Slimane -Moron  Path 

The  results  of  transmission  loss  measurements  over  this 

8  9  / 

path  have  been  described  previously— 1 — '.  In  order  to  improve  the 
prediction  for  the  operation  in  the  700-900  Me  range  used  as  examples 
in  Ref.  1,  the  17  days  of  measurements  over  exactly  the  same  path 
(Sidi  Slimane  Relay  Center  to  Moron  Site  A)  on  1046  Me  were  utilized. 
For  the  purpose  of  this  example  it  is  assumed  that  the  difference  in 
frequency  is  not  great  enough  to  make  the  paths  significantly  different, 
i.e.  ,  pak  is  assumed  to  be  equal  to  1.  The  measured  distribution  of 
hourly  median  values  of  basic  transmission  loss  is  shown  by  the  dotted 
line  on  Fig.  7,  designated  as  The  measurements  were  made  in 

February  of  1957  during  all  hours  and  thus  represent  time  block  group 
F  (Time  Blocks  1,  2,  3  and  8). 

In  accordance  with  (16)  of  Sec.  3,  above,  the  measured 
values  are  adjusted  to  all  day-all  year  values  by  subtracting  the  term 
AV^p(p,  9).  This  term,  given  by  (15),  is  the  difference  between  the 
V(p,  0)  values  for  all  day -all  year  (read  from  Fig.  4  at  the  path 

angular  distance  value  0  =  40.  3  milliradians) ,  and  the  V(p,  0)  values 

8,9/  A.  P.  Barsis  and  R.  S.  Kirby,  "Long  distance  tropospheric 
propagation  measurements  between  North  Africa  and  Spain,  " 

NBS  Reports  5067  and  5524,  May  15  and  October  15,  1957 
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for  all  hours -winter  (read  from  Fig.  2  at  the  same  0  =  40.  3  mill!  - 
radians).  The  adjusted  values,  Lbo(p),  are  shown  by  the  broken  line 
on  Fig.  7.  This  curve  represents  the  observed  values.  The  weight 
assigned  to  the  median  L^Q(50)  is  obtained  from  (9),  using  Figs.  5  and 
6.  For  n  =  17  days,  cro^(50)  =  27.  1  (db)^,  and  for  0  =  40.  3  milliradians, 
f ( 0)  =  0.  398.  Thus  <r^(50)  is  determined  to  be  10.  79  (db)^.  Values  of 
°"q{p)  are  now  determined  for  other  percentage  values  in  accordance 
with  (17)  using  the  differences  of  the  appropriate  V(p,0)  values  from 
Figs.  2  and  4,  as  before.  To  each  cr^(p)  a  term  cr^  =4.0  (db)  is 
added  in  order  to  allow  for  equipment  (reading  and  calibration)  errors. 
The  resultant  values  crQe(p)  represent  the  uncertainty  assigned  to  each 
observed  value  L^^p),  as  adjusted  to  all  day-all  year  values.  The 
weight  wQ(p)  assigned  to  each  value  Lbo(p)  is  the  reciprocal  of  0"Og(p). 

For  the  same  1046  Me  frequency,  the  Time  Block  2  median 
(of  hourly  medians)  basic  transmission  loss  value  is  read  from  Fig.  2 
of  Ref.  1,  as  the  sample  calculations  of  Ref.  1  refer  to  the  same  path. 
This  value,  Lbm  =  214.  7  db,  together  with  the  V(p,  0)  values  of  Fig.  4, 
produces  the  calculated  distribution  of  hourly  medians  for  all  day -all 
year; 

Lbc(p,  0)  =  Lbm  -  V(p,  0)  (22) 

The  calculated  values  Lbc(p,  0)  are  shown  by  the  dash-dotted 
line  on  Fig.  7.  They  may  be  redesignated  as  L-bc(P)>  as  the  angular 
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distance  9  has  now  been  fixed.  The  prediction  uncertainty  crc(p) 
associated  with  each  L^c(p)  is  given  by  curve  (a)  of  Fig.  1,  as  there 
is  no  significant  contribution  of  the  diffracted  field  over  this  path,  and, 
prior  to  combining  predictions  and  measurements,  the  system  error 
(rr  (representing  principally  the  uncertainties  in  estimating  the  system 
parameters  F  and  R)  does  not  enter.  The  weights  wc(p)  =  l/cr^p)  are 
shown  on  Fig.  1 . 

Now  the  weighted  average  of  L^c(p)  and  L|=)0(p)  may  be 
obtained  in  accordance  with  (la),  resulting  in  the  distribution  of  hourly 
medians  L^^.(p),  the  improved  prediction.  This  is  the  solid  curve  on 
Fig.  7,  and  the  corresponding  weights  w^(p)  are  obtained  by  adding 
wc(p)  and  wQ(p). 

The  performance  of  a  communication  system  utilizing  this 
path  may  now  be  predicted  on  the  basis  of  three  methods:  first,  using 
the  calculated  values  alone;  second,  using  only  observed  values;  and 
third,  using  the  weighted  average  of  the  two,  which  constitutes  the 
improved  prediction. 

The  first  case  has  already  been  discussed  in  Ref.  1.  For 

this  case  the  distribution  L^c(p)  of  Fig.  7  is  used  with  the  prediction 

uncertainty  given  by  the  dash-dotted  o'rc(p)  curve  of  Fig.  8,  which  is 

the  same  curve  as  the  (b)  curve  of  Fig.  1.  The  system  error  crr  is 

included  in  tr  (p). 

rc 
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crv|(p)  =  <d(p)  +  ^  (23) 

For  the  second  case  the  observed  distribution  (adjusted  to 

all  day -all  year)  L^Q(p)  of  Fi.g.  7  is  used,  together  with  the  prediction 

uncertainty  o'Qr(p)  shown  by  the  dashed  curve  on  Fig.  3.  As  the 

prediction  applies  to  a  complete  system,  cror(p)  also  contains  the 

system  error  a  : 

r 

^or^P)  =  °oe(p)  +  ar  (23a) 

Finally,  the  improved  prediction  curve  L^t(p)  of  Fig.  7  is 
used  together  with  the  prediction  uncertainty  o’tr(p)  shown  by  the  solid 
curve  on  Fig.  8.  Again,  a^r(p)  contains  the  system  error  crr,  and  is 
given  by 

°-tr(P>  =  <rf  (P)  +  \  (23b) 

2 

where  cr^  (p)  is  the  reciprocal  of  the  sum  of  the  weights  wQ(p)  and 
we(p)  in  accordance  with  (20). 

Fig.  8  shows  that  the  prediction  uncertainty  is  appreciably 
reduced  by  the  use  of  these  measurements.  Note  that  the  distribution 
of  the  observed  hourly  median  values  shown  on  Fig.  7  indicates  sub¬ 
stantially  lower  fields  at  the  high  percentage  values  than  would  be 
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expected  from  the  values  calculated  on  the  basis  of  the  V(p,  9)  curves 
of  Fig.  4.  * 

As  the  last  step,  service  probability  curves  will  be  cal¬ 
culated  in  accordance  with  the  methods  of  Ref.  1.  In  this  particular 
case,  the  effects  of  higher  transmission  loss  and  lower  prediction  un¬ 
certainty  values  tend  to  balance  each  other,  as  seen  by  the  final 
comparison  of  service  probabilities  for  99  or  99.  9%  of  all  hours  for 
the  three  cases. 

It  was  explained  in  Ref.  1  that  the  system  power  P[F(t),p] 
required  to  provide  the  specified  grade  of  service  or  better  during  p % 
of  all  hours  with  the  probability  F(t)  is  related  to  the  expected  value  of 
system  power  P(0.  5,  p)  by 

P[F(t),p]  =  P(0.  5,p)  +  t  .  cr(p)  (24) 

where  P(0.  5,  p)  here  relates  to  the  three  curves  of  basic  transmission 
loss  values  L]3c(p),  ^  (p),  and  L^t(p),  and  cr(p)  is  identical  to  either 
crrr(P)«  o'  (p),  or  a,  (p)  depending  on  the  three  cases  considered. 

For  the  same  system  as  used  in  Ref.  1  (60  ft.  parabolic  reflectors, 
quadruple  diversity  FM,  850  Me  carrier  frequency),  and  the  same 

specified  grade  of  service  (24  voice  channels,  several  of  which  are 

*  It  should  be  noted  that  the  original  measured  distribution  (Fig.  3 
of  Ref.  7)  has  been  smoothed  and  extrapolated  in  order  to  facilitate 
the  calculations. 
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loaded  with  standard  frequency -shift  teletype  signals  with  less  than 
1  in  10, 000  character  errors)  it  can  be  shown  by  substitution  into 
(19)  that  the  required  transmitter  power  in  decibels  above  one  watt 
is  given  by  the  hourly  median  basic  transmission  loss  minus 
197  db.  Thus  a  fixed  transmitter  power  Pq  -  1  kw  (30  dbw)  corresponds 
to  227  db  basic  transmission  loss  on  the  curves  of  Fig.  7,  and  PQ  =  10  kw 
(40  dbw)  corresponds  to  =  237  db.  Substituting  PQ  for  P[  F(t),  p]  in 
(23)  and  solving  for  t, 


t  = 


Po  -  P(0.  5,  p) 


o-(P) 


(25) 


and  from  t,  the  service  probability  F(t)  may  be  found  from  appropriate 
probability  tables  10/. 

For  each  of  the  three  cases,  F(t)  may  be  now  determined  as 
a  function  of  p,  with  P(0.  5,  p)  and  or(p)  read  from  the  appropriate  curve 
of  Figs.  7  and  8,  and  for  Pq  fixed  at  1  or  10  kw.  The  resultant  curves 
of  service  probability  are  shown  on  Fig.  9. 

The  curves  on  Fig.  9  show  the  anticipated  results,  namely 
that  the  decrease  in  the  prediction  uncertainty  is  not  sufficient  to 
compensate  for  the  effect  of  the  low  values  of  field  measured  during  an 

appreciable  percentage  of  the  time.  Only  on  the  basis  of  a  0.  99 

10/C.  A.  Bennett  and  N.  L.  Franklin,  "Statistical  analysis  in 

chemistry  and  the  chemical  industry,  "  John  Wiley  and  Sons,  Inc.  , 
New  York,  1954.  See  p.  8  6  and  Table  I  of  the  Appendix 
(pp.  689-693) 
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service  probability  (only  1  out  of  100  chances  of  failure)  do  the  mea¬ 
surements  indicate  a  greater  time  availability  of  service  than  is  in¬ 
dicated  by  the  prediction  alone.  This  applies  to  the  use  of  1  kw  trans¬ 
mitter  power.  The  analysis  suggests  that  for  this  particular  example 
the  use  of  the  calculated  values  alone  would  probably  lead  to  results 
that  were  too  optimistic,  especially  if  lower  values  of  service  probabi¬ 
lity  (0.  5  to  0.  8)  had  been  the  basis  for  the  estimates. 

A  more  direct  comparison  may  be  made  by  fixing  F(t)  and  p, 
and  solving  for  PQ  in  (25);  this  shows  the  transmitter  power  necessary 
to  provide  the  specified  grade  of  service,  or  better,  during  p%  of  all 
hours  with  a  given  probability  of  success  F(t).  If  this  is  done  for  the 
sample  path  and  the  system  considered,  the  following  results  are 
obtained: 

Table  IV 


Transmitter  Power  in  Decibels 
Relative  to  1  Watt 


(a)  F(t)  =  0.  95,  p  =  99% 

Calculated  only  38.  5 

Observed  only  45.  1 

Improved  prediction  39.  7 

(b)  F(t)  =  0.  99,  p  =  99% 

Calculated  only  40.  3 

Observed  only  48.  9 

Improved  prediction  42.  9 

(c)  F(t)  =  0.  99,  p  =  99.  9% 

Calculated  only  48.  1 

Observed  only  55.  1 

Improved  prediction  48.  2 
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In  Table  IV,  (a)  corresponds  to  a  0.95  probability  (1  out  of  20  chances 
of  failure)  to  obtain  the  specified  grade  of  service  or  better  during 
99%  of  all  hours.  For  (b)  the  probability  is  increased  to  0.99  (only 
1  out  of  100  chances  of  failure),  and  for  (c)  the  percentage  of  hours  is 
increased  to  99.9%,  representing  together  with  the  0.99  service 
probability  a  virtually  perfect  circuit.  In  each  case  the  required  trans 
mitter  power  would  be  substantially  higher  if  based  on  the  analysis  of 
the  limited  observations  only.  The  results,  based  on  the  predicted 
values  improved  by  the  measurement,  are  not  very  different  from  the 
results  based  on  the  use  of  the  calculated  predicted  values  alone. 

For  the  path  and  the  system  under  discussion,  therefore, 
the  measurements  did  not  add  significantly  to  the  information  gained 
from  theoretical  considerations  alone,  and  estimates  based  on  the 
use  of  the  measurement  alone  would  have  been  entirely  too  pessimistic 
These  conclusions,  of  course,  apply  to  service  fields,  where  a  high 
degree  of  reliability  during  the  greatest  percentage  of  hours  of  the 
year  is  desired.  The  findings  would  have  been  quite  different  in  the 
consideration  of  nuisance  or  interference  fields  for  this  path. 

(b)  Prediction  Improvement  for  the  St.  Louis -Urbana  Path 

Field  strength  measurements  on  broadcasting  station 
KXOK-FM  were  made  over  the  period  July  1950  to  June  1952  at 
Urbana,  Illinois.  This  is  one  of  the  paths  listed  in  Appendix  IV  of  a 
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previous  comprehensive  paper  on  the  use  of  angular  distance  in  trans¬ 
mission  loss  calculations^/;  the  path  distance  is  146.  5  miles,  the 
angular  distance  19.3  milliradians ,  and  the  operating' frequency 
93.  7  Me. 

In  order  to  investigate  prediction  improvement  as  a  function 
of  the  length  of  the  measurement  period,  the  following  nominal  time 
periods  were  selected  from  the  KXOK  data: 


Nov.  1  -  7,  1950  .  7  days 

Nov.  1  -  15,  1950  .  1  5  days 

Nov.  1  -  30,  1950  .  30  days 

Nov.  1  -  Dec.  31,  1950  .  60  days 

Nov.  1,  1950  -  Apr.  30,  1951  .  180  days 

Nov.  1,  1950  -  Oct.  31,  1951  .  365  days 


For  each  period  all  available  hourly  medians  were  used,  thus 
all  but  the  365-day  period  represent  the  winter  time  block  group  F,  and 
the  365-day  period  represents  all  hours  of  the  year.  It  may  be  assumed 
that  a  measurement  program  was  started  on  Nov.  1,  1950,  and 
successively  longer  periods  of  measurement  were  used  for  the  purpose 
of  prediction  improvement.  Fig.  10  shows  the  measured  distributions 
of  hourly  medians  of  basic  transmission  loss  for  the  first  four  periods 
tabulated  above,  and  Fig.  11  the  measured  and  theoretical  distribu¬ 
tions  for  the  full  time  block  group  F  (all  hours,  winter),  and  for  all 
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hours  of  the  year.  The  theoretical  distributions  were  obtained  from 
Figs.  2  and  4  using  the  calculated  Time  Block  2  median  basic  trans¬ 
mission  loss  value  of  179.2  db  as  the  reference  value.  L,  .  This,  in 

bm 

turn,  was  calculated  in  accordance  with  the  method  given  in  Ref.  5. 

Using  the  same  procedure  as  in  the  previous  example,  mea¬ 
sured  distributions  were  adjusted  to  all  day-all  year  values,  and  com¬ 
bined  with  the  calculated  all  day -all  year  distribution. 

The  following  table  lists  0^(50)  versus  n  derived  from 
Figs.  5  and  6  for  this  path,  with  f(9)  -  0.  93: 

Table  V 


Number  of 
days 

cr  2(50) 
oo'  ' 

^(50) 

n 

db 

db 

7 

38.  85 

36.  1 

15 

28.  7 

26.  7 

30 

21.2 

19.  7 

60 

15.  7 

14.  6 

180 

9.4 

8.  75 

365 

6.  7 

6.  25 

The  resultant  improved  prediction  curves  of  the  expected 
distributions  of  hourly  medians  are  shown  on  Fig.  12  taking  into 
account  all  six  measurement  periods.  The  appropriate  standard  devia¬ 
tions  (r^r  characterizing  the  remaining  prediction  uncertainties  are 
shown  on  Fig.  13.  Each  of  the  curves  is  labeled  with  the  number  of 
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days  contained  in  the  measurement  period  utilized.  The  standard 
deviation  in  each  case  contains  the  component  crr>  as  in  the 
previous  example,  in  order  to  take  system  performance  uncertainties 
into  account.  It  is  to  be  kept  in  mind  that  the  measurements  would 
ultimately  be  used  to  estimate  the  performance  of  a  communication 
system. 

For  the  purpose  of  illustrating  the  problem  it  is  sufficient 
to  assume  that  one  single  value  of  hourly  median  basic  transmission 
loss  on  the  curves  of  Figs.  1  1  or  1Z  characterizes  a  system  with  a 
fixed  transmitter  power  value  P^.  The  value  chosen  for  the  example 
is  190  db,  and  it  is  taken  as  the  basis  for  service  probability  calcula¬ 
tions.  For  each  of  the  L^^p)  curves  of  Fig.  12,  the  appropriate 
o"tr(p)  curve  (for  the  same  number  of  days)  of  Fig.  13  is  used.  The 
resulting  service  probability  curves  are  shown  on  Fig.  14;  they  are 
again  labelled  with  the  number  of  days  to  which  they  apply.  Fig.  14 
also  includes  the  service  probability  curve  for  calculated  values 
alone  (without  taking  into  account  the  measurements).  These  are 
based  on  the  same  190  db  value  of  basic  transmission  loss  characterizing 
the  assumed  system,  and  were  developed  using  the  calculated  all  day- 
all  year  distributions  of  Fig.  11  with  the  a^.c(p)  standard  deviation  curve 
of  Fig.  8. 
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As  in  the  previous  example,  any  of  the  service  probability 
curves  in  Fig.  14  is  affected  not  only  by  the  remaining  prediction  un¬ 
certainty  expressed  by  the  standard  deviation  ertr(p),  but  also  by  the 
shape  of  the  appropriate  L^^p)  curve,  and  thus  by  the  measured 
distribution  of  hourly  medians.  Figs.  10,  11  and  12  show  that  the 
first  60  days  of  measurement  did  not  produce  enough  high  trans¬ 
mission  loss  (low  field  strength)  values  to  pull  the  curves  down  suffi¬ 
ciently  in  the  high  percentage  range.  Only  for  the  longer  measure¬ 
ment  periods  the  improved  prediction  takes  into  account  those  low 
field  strength  values.  The  standard  deviation  (Tj-r(p),  of  course, 
decreases  systematically  inversely  to  the  length  of  the  measurement 
period,  but  does  not  always  compensate  for  the  low  field  strength 
values  observed  during  the  extended  measurement  period. 

Summarizing  some  of  these  observations,  two  things  happen 
to  a  service  probability  curve  when  additional  measurements  are 
taken  into  account.  The  slope  of  the  curve  decreases  with  the  variance 
(^(p),  as  more  and  more  observations  are  taken  into  account.  There  is 
also  a  change  in  level  of  the  service  probability  curves,  depending  on 
the  behavior  of  the  observed  data.  This  change  in  level  may  be  either 
up  or  down  if  additional  data  are  added,  as  illustrated  by  the  KXOK 
example.  The  combined  effect  of  a  change  in  slope  and  a  change  in 
level  may  be  either  an  increase  or  a  decrease  in  service  probability 
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at  a  given  value  of  time  availability.  The  range  of  oscillation  of  these 
estimates  of  service  probability  will  decrease  as  the  number  of  mea¬ 
surements  is  increased.  The  Sidi  Slimane -Moron  example  may  very 
well  represent  an  extreme  departure  from  this  expected  trend.  If,  in 
that  case,  additional  measurements  had  been  performed  comparable 
in  extent  to  the  KXOK  example,  a  similar  pattern  for  prediction  im¬ 
provement  as  a  function  of  the  length  of  the  measurement  period  might 
have  been  recognized. 

Returning  to  the  KXOK  example  and  Fig.  14,  it  is  seen  that 
for  the  service  probability  value  of  0.  95  any  of  the  curves  based  on  the 
improved  prediction  indicate  a  higher  value  of  time  availability  than  the 
curve  based  on  calculated  values  alone. 

As  in  the  previous  example,  a  comparison  of  basic  transmis¬ 
sion  loss  values  (which  directly  correspond  to  system  power)  may  be 
made  for  various  lengths  of  the  measurement  period  at  fixed  values  of 
p  and  F(t).  The  results  are  outlined  in  the  following  table  for  F(t)  =  0.  95 
and  0.  99,  and  a  time  availability  p  of  99%  and  99.  9%  of  all  hours. 
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Table  VI 


Equivalent  Basic  Transmission  Loss  in  Decibels 


P  =  99% 

F (t)  =  0.  95 

P  =  99% 
F(t)  =  0.  99 

P  =  99.  9% 
F(t)  =  0.  99 

Based  on  calculated 

values  only 

Based  on  improved 

predic  tion 

207.  3 

212.  2 

220.  6 

by  measurement. 

201.  3 

205.  4 

211.  7 

14  days 

200.  2 

204.  1 

211.0 

1  month 

200.  5 

204.  3 

211.5 

Z  months 

200.  0 

203.  7 

210.  8 

6  months 

200.  7 

204.  4 

211.8 

1  year 

200.  9 

204.  5 

211.  8 

This  particular  example  shows  that  the  method  of  combining 
calculated  and  observed  values  may  provide  a  substantial  change  in 
the  prediction,  with  a  6  to  8  db  reduction  in  transmitter  power  required 
to  provide  the  specified  grade  of  service  during  almost  all  hours  of 
the  year  with  a  very  high  probability  of  success.  In  this  case,  it  also 
appears  that  only  7  days  of  measurements  are  adequate  to  give  this 
result,  and  the  extension  of  measurement  over  the  period  of  one  full 
year  does  not  provide  much  further  change.  On  the  other  hand,  some 
other  seven-day  period  could  have  been  selected  with  different  results. 

Of  course,  it  must  be  remembered  that  all  the  measured 
field  strength  values  for  90%  of  all  hours  or  more  were  higher  than 
the  calculated  ones.  This  is  the  opposite  condition  in  comparison  to 
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the  other  example.  It  may  be  concluded,  generalizing  from  this 
example  only,  that  a  good  indication  of  the  prediction  improvement  may 
be  obtained  from  only  seven  days  of  observation,  while  a  larger  mea¬ 
surement  period  serves  only  to  confirm  the  result. 

6.  General  Conclusions 

As  an  extension  of  a  method  permitting  evaluation  of  the  probabil¬ 
ity  of  success  of  a  tropospheric  communication  system  on  the  basis  of 
empirical  parameters  obtained  from  an  extensive  and  very  general 
measurement  program,  it  has  been  shown  how  this  prediction  method 
may  be  modified  to  take  into  account  transmission  loss  measurements 
over  exactly  the  same  path  or  a  similar  path  for  which  the  final  system 
is  intended.  Of  the  two  examples  used  to  demonstrate  the  procedure, 
one  indicates  a  substantial  increase  in  time  availability  at  a  given  level 
of  prediction  uncertainty  as  estimated  using  only  a  short  period  of 
measurements,  while  the  other  one  does  not  indicate  any  significant 
change,  taking  account  of  measurements.  In  both  examples  it  is  shown 
that  one  has  to  be  very  careful  in  drawing  conclusions  from  the  results 
of  measurement;  the  reduction  in  the  prediction  uncertainty  does  not 
always  compensate  for  extreme  departures  of  observations  from  the 
calculated  values.  The  "improved  prediction"  obtained  from  forming  a 
weighted  average  of  calculated  and  observed  values  will  always  result 
in  a  reduction  of  the  prediction  uncertainty,  because  the  resultant 
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weight  is  greater  than  each  of  the  individual  weights  assigned  to  either 
the  calculated  or  the  observed  values.  However,  the  actual  departure 
of  calculated  from  observed  field  strength  (or  transmission  loss  values) 
plays  an  important  part,  as  shown  by  the  examples. 

In  some  cases  the  designer  or  engineer  wants  to  know  whether 

to  place  reliance  on  an  available  set  of  measurements,  or  on  a 

prediction  formula.  Thus  it  may  be  of  interest  to  determine  the 

expected  number  of  days  of  measurement  necessary  to  result  in  a 

prediction  uncertainty  equal  to  the  prediction  uncertainty  associated 

with  the  calculated  value.  In  terms  of  the  preceding  discussions,  this 

2 

means  setting  (17)  equal  to  <r  (p)  as  given  on  Fig.  1  and  solving  for  n, 
the  number  of  days  of  measurements,  as  a  function  of  the  angular 
distance  8,  and  the  time  availability  p  in  percent  of  all  hours. 

Several  factors  have  to  be  recognized  immediately: 

1.  A  certain  minimum  number  of  hours  (or  days)  is  necessary 
to  define  a  given  percentage  value  for  a  distribution  of  hourly  medians. 
As  an  example,  it  is  necessary  to  have  at  least  50  hourly  medians  in 
order  to  define  a  99%  value. 

2.  If  not  more  than  180  consecutive  days  of  measurements  are 
available,  they  may  fall  entirely  within  time  block  group  F  (winter)  or 
S  (summer).  In  this  case,  allowance  has  to  be  made  in  order  to  apply 
the  results  to  all  day -all  year  conditions. 
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3.  The  variance  of  the  observed  hourly  medians  is  strongly 
dependent  on  the  angular  distance,  9,  whereas  the  variance  of  the 
calculated  hourly  medians  has  so  far  been  found  to  be  practically 
independent  of  0,  at  least  for  9>  10  milliradians,  where  the  contribution 
of  the  diffracted  field  is  negligible. 

Thus  we  may  write  in  accordance  with  (17),  averaging  the  effect 
of  the  winter  and  the  summer  time  block  groups. 

0-o|(P)  =  <r2(p)  +  <rf  =  <ro|(50)  f(0)  +  0.  09[  (AVA)2 


+  1/2(AVaf)2  +  1/2(AVas)2]  (26) 

where: 

AVA  =  VA(P’e)’  '  VA<50'8>  <27> 

AVaf  =  VA(p,  6)  -  VF(p,  0)  (27a) 

AVas  =  VA(p.  0)  -  Vs(p,  0)  (27b) 


The  subscripts  A,  F,  and  S  refer  to  all  day-all  year,  the  winter  time 
block  group  and  the  summer  time  block  group,  respectively,  with  the 
V(p,  0)  values  read  from  Figs.  2,  3  or  4. 

(26)  is  solved  for  0"Oo(5O)  with  p  =  50,  90,  99  and  99.  9%  as  a 

2  2  2 
function  of  9,  and  cr  assumed  to  be  4  db  ,  by  setting  o-  P(p)  equal  to 

°"q(P)»  the  latter  being  read  from  Fig.  1  for  the  various  values  of  p 

used.  The  number  of  days,  n,  is  then  obtained  from  Fig.  5  as  a 
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function  of  cr  ^(50). 

oo 

Results  are  plotted  on  Fig.  15,  with  the  number  of  days,  n,  as 
ordinate  versus  the  angular  distance  9  plotted  as  the  abscissa.  The 
curves  shown  are  labelled  with  time  availability  in  per  cent  of  all  hours. 
The  results  are  interpreted  in  the  following  way;  taking  into  account 
equipment  (calibration  and  reading)  errors,  it  takes,  as  an  example, 

10  days  of  continuous  measurements  to  obtain  a  value  of  prediction 
uncertainty  cr0e(99%)  at  9  =  29  milliradians  for  99%  of  all  hours,  which 
is  equal  to  the  prediction  uncertainty  (r^(99%)  assigned  to  the  calculated 
value.  In  this  particular  example,  the  <r  will  be  smaller  than  the 
croe(99%)  for  any  shorter  measurement  period;  thus  the  calculated 
distribution  of  hourly  medians  would  constitute  a  better  estimate  than  a 
distribution  obtained  from  observations  over  a  shorter  time  period 
than  19  days.  If  the  measurement  period  is  longer  than  19  days,  the 
opposite  is  true,  but  in  any  event  a  still  better  estimate  is  provided  by 
combining  calculated  and  observed  values  in  the  manner  described  in 
this  paper. 

Fig.  15  should  not  be  interpreted  as  indicating  where  and  for  how 
long  measurements  should  be  made.  This  question  still  depends  on  the 
relative  cost  of  measurements  versus,  say,  the  cost  of  more  trans¬ 
mitter  power  or  larger  antennas.  If,  for  instance,  the  calculations  for 
a  particular  path  result  in  the  requirement  of  19  kw  transmitters  and 
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60  ft.  parabolic  reflectors,  which  would  provide  the  specified  grade  of 
service  or  better  with  a  probability  of  0.  95  during  99%  of  all  hours,  and 
if  two  weeks  of  measurements  might  be  expected  to  show  that  28  ft. 
antennas  would  provide  a  similar  service  probability,  the  designer  will 
have  to  balance  the  cost  of  the  measurements  against  the  cost  of  in¬ 
stalling  and  maintaining  the  larger  antennas  with  due  regard  to  the  fact 
that  the  larger  antennas  would,  in  any  case,  provide  a  better  service. 

Of  course,  there  is  always  a  possibility  that  measurements  will  show 
that  more  power,  or  larger  antennas  are  required.  In  general,  it  is 
quite  likely  that  measurements  over  a  particular  path  will  only  serve 
to  confirm  the  estimates  made  solely  on  the  basis  of  calculations. 

Although  the  basic  concept  of  service  probability  and  the  methods 
for  prediction  and  prediction  improvement  will  always  be  valid,  some 
changes  in  the  detailed  procedures  are  expected  on  the  basis  of  addi¬ 
tional  data  evaluation  and  interpretation.  Work  is  in  progress  now  to 
re-evaluate  the  long-term  variability  of  hourly  medians,  and  improve 
their  reliability  by  allowing  for  the  effect  of  frequency,  especially  at 
short  distances.  Also,  it  may  turn  out  to  be  slightly  better  to  replace 
the  dependence  on  the  angular  distance  9  by  a  dependence  on  antenna 
height  and  linear  path  distance.  A  better  understanding  of  the  factors 
underlying  the  long-term  variability  will  ultimately  result  in  a  signi¬ 
ficant  reduction  of  the  prediction  uncertainty. 
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APPENDIX  I 
Definitions  of  Terms 


This  table  contains  most  of  the  symbols  and  terms  used  in  Secs.  2, 

3,  and  4  of  this  paper. 

Lbc(P)  -  calculated  hourly  median  basic  transmission  loss  not 
exceeded*  during  at  least  p %  of  all  hours  within  the 
time  block,  or  the  time  block  combinations  con¬ 
sidered.  If  desired,  p  may  have  a  subscript  de¬ 
noting  the  number  of  the  time  block,  and  subscripts 
F  and  S  may  be  used  to  denote  combined  Time 
Blocks  1,  2,  3,  8  (all  hours,  winter,  and  4,  5,  6, 

7  (all  hours,  summer),  respectively.  Thus,  L^m 

from  Sec.  3  of  Ref.  1  may  be  designated 

For  the  definition  of  the  time  blocks,  see  p.  11. 


°c(p) 
wc(p)  = 


standard  deviation  characterizing  the  prediction  un¬ 
certainty  of  a  calculated  value  L^^p) 

-  -  weight  assigned  to  the  calculated  value  L^c(p). 


(p)  -  observed  hourly  median  basic  transmission  loss 

value  not  exceeded  during  at  least  p%  of  all  hours 
*In  terms  of  field  strength,  L^c(p)  corresponds  to  the  hourly  median 
field  strength  value  exceeded  during  at  least  p %  of  all  hours. 
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within  a  time  period  which  is  defined  in  the  same 


way  as  explained  for  L^ip)  above. 


standard  deviation  characterizing  the  uncertainty 


errors  introduced  by  the  measuring  equipment  and 
procedures. 


standard  deviation  characterizing  errors  introduced 


a 


e 


by  measuring  equipment  and  procedures. 


wQ(p) 


-  weight  assigned  to  the  observed  value  Lbo<P>- 


taking  into  account  c : 


standard  deviation  characterizing  the  uncertainty  in 


a 


r 


estimating  the  system  parameters  F  (noise  figure) 
and  R  (median  r.m.  s.  carrier  to  r.m.  s.  noise  ratio). 


<rAV(P)  =  k(AV)  '  standard  deviation  associated  with  the  distri¬ 


bution  of  differences  in  the  long  term  variable, 
V(p,  0);  taken  for  a  given  value  of  0,  but  for 
various  percentage  levels  p.  The  values  of 
V(p,  0)  which  define  AV  may  be  taken  in  the  same, 
or  in  different  time  blocks,  or  time  block  groups. 
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aor(P)  "  standard  deviation  characterizing  the  prediction 
uncertainty  based  on  measurements  only,  but  in¬ 
cluding  a  (see  above)  in  order  to  be  applicable  to 
the  evaluation  of  systems 


(P)  + 


Lbt(F) 


°-t(p) 


weighted  average  of  predicted  and  observed  trans¬ 
mission  loss  values  Lbc(P)  and  Lbo(p)  (see  above), 
standard  deviation  characterizing  the  remaining 
prediction  uncertainty  associated  with  the  weighted 
average  Lbt(p) 

— \ =  wt(p)  =  wQ(p)  +  wc(p) 

<rt2(P) 


O' t^ (p)  -  standard  deviation  characterizing  the  remaining 

prediction  uncertainty  of  Lbt(p)  as  based  on  both 
prediction  and  measurements,  and  including  o' 

(see  above)  in  order  to  be  applicable  to  the  evalua¬ 
tion  of  systems 


°-tr(P)  =  °-t(P)  +  o-r 
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APPENDIX  II 

Estimating  the  Prediction  Uncertainty  on  the  Basis  of 

Measurements  Over  a  Different  Path 


Starting  from  (2)  in  Sec.  2  of  the  report: 


^tb  ^ta  ^  ^cb  ~  ^ca^ 


(II-l) 


This  means  that  the  improved  prediction  for  Path  b  may  be  ex¬ 
pressed  as  the  improved  prediction  for  Path  a  plus  the  difference  in 
the  calculated  values  for  the  two  paths.  The  improved  prediction  Lta 
for  Path  a  was  given  by  (1)  of  Sec.  2: 


Jta 


woa  Loa  +  wCa  Lca 

wo  +  wca 


(II-2) 


7 

Its  uncertainty  may  be  expressed  by  the  variance  or^,  which  is  the 
reciprocal  of  the  sum  of  the  weights  wQa  +  wca  (denominator  of  II — 2 ) , 
The  predicted  values  L  _  and  L  .  include  the  true  values  L_ 

Cq  CD  B. 

and  plus  the  prediction  errors  ea  and 


L  -  L  +  c  ;  L  ,  -  Lu  +  e , 
ca  a  a  cb  13  b 


(H-3) 


The  observed  value  Loa  appearing  in  (II - 2)  is  the  sum  of  the 


true  value  L  and  a  measurement  error  6_ 

3,  cL 
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L 

oa 


=  L  +6 

a.  a 


(II  -4) 


The  true  values  La  and  are  constants,  and  the  errors  &a,  c  , 
and  e  ^  are  statistical  variables  with  means  of  zero  for  each,  and  their 
variances  given  by  cr  ^  =  l/w-a,  cr  2  =  l/\v  .  and  cr  ?  =  l/w  ,  , 
re  spectively. 

Inasmuch  as  the  expected  values  of  L^a  and  L^,  as  well  as  Lca 
and  are  La  and  ,  (II- 1)  reduces  to  the  identity 

Lb  -  La  =  Lb  -  La  <n-5» 


when  written  in  terms  of  expected  values. 

The  variance  of  the  linear  combination  of  the  three  terms  in  (II-l) 
is  given  byi/ 

2  ?  7  2 

crfi  -  c.r'  +  cr  “  +  cr  ,  -  2c r,  cr  p.  _  +  2<r,  cr  ,  p ,  „ 

tt>  ta  ca  cb  ta  ca  K1Z  ta  cb 


-  2cr  cr  ,  p„  _ 
ca  cb  *23 


(II-6) 


In  (II -6)  above  p  designates  the  correlation  coefficient  between  the 
improved  prediction  L  and  the  calculated  value  L  for  Path  a,  p.  _ 
the  correlation  coefficient  between  the  improved  prediction  L.  for 

Z3i 

Path  a  and  the  calculated  value  L  ,  for  Path  b,  and  p_,0  the  correlation 

c  b  rZ3 

coefficient  between  the  calculated  value  L  for  Path  a,  and  the  cal- 

ca 

1/  See,  for  instance,  C.  A.  Bennett  and  N.  L.  Franklin,  "Statistical 
analysis  in  chemistry  and  the  chemical  industry,  "  John  Wiley  and 
Sons,  Inc.  ,  New  York,  1954.  See  p.  50. 


-50- 


culated  value  L  ,  over  Path  b. 
cb 

From  Ref.  1  (see  p.  37)  the  correlation  coefficient  px^.  between 
the  variables  x  and  y  may  be  written  as 

„  _  E(xy)  -  E(x)  E(y) 


xy 


tr  cr 
x 


(H-7) 


y 


where  E(x),  E(y),  and  E(xy)  denote  here  expected  values  of  the 
variables  x,  y,  and  of  their  product  xy. 

Thus , 

_  E<LtaLca>  "  E<Lta>  E<Lca> 


12 


(ir-8) 


O'.--,  O’ 
ta  c  a 


Substituting  (II-2)  and  ( II - 3)  into  ( II -8) ,  and  noting  that 


w  +  w 
oa  oc 


crta,  we  obtain 


12  cr.  tr  „ 
ta  ca 


E 


2  ? 
cr.  (w__  L_  _  L  +  w  L  „) 

ta  '  oa  oa  ca  ca  ca 


-  E 


a./  (w  L  +  w  L  ) 
ta  oa  0a  ca  ca' 


“ 

** 

cr,  ^  (w  L  _  +  w  L  ) 

ta  '  oa  oa  ca  ca' 

.  E 

L 

ca_ 

J 

(II -9) 
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Expanding  the  various  terms  separately, 


Loa  LCa 


=  E 


<La  -  6a>  'La  "  €  a) 


=  E 


Lj  '  La  “  e 

a  a  a  a  a  a-  a 


=  L 


a  * 


(II- 10) 


As  E(6a€a)  equals  zero  for  the  reasonable  assumption  that  the  pre¬ 
diction  and  measurement  errors  are  not  correlated, 


E 

[L ; 

=  E 

L2  -  2L„  €  -f  e  2 

a  a  a  a 

-  T  2  4-  n-  2 

-  La  +  aca 


(II-H) 


<r.  f(w  +  w  L  ) 

tav  oa  oa  ca  ca 


L 


ca 


=  La 


(11-12) 


Finally, 


p12 


ata  ^ca 


ta 


w  L2  +  w  (L^  +  cr  2) 
oa  a  ca  '  a  ca' 


-  L1 


cr  <r 
ta  ca 


1  J  'tiKa  +  wca)  +  -ta  ‘  La}  =  ^  I11’13) 


Similarly, 

E(LtaLcb)  -  E(Lta)  •  E(Lcb) 

^l2  ^ta  ^cb 


(11-14) 
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Pl3  'rta°cb 


K 


^ta  ^woa  ^oa  ^cb  +  wca  ^ca 


-  E 


ataKaLoa  +  wcaLca> 


L 


cb 


]} 


(11-15) 


Analogous  to  (II  -10) , 


Loa  -^cb 


=  E 


<La  -  5a)  (Lb  -  €b) 


=  LaLb,  (II-16) 


and 


Lc a  Lc^  -  E 


(La  “  €  a)  (^b  e 


-  La  Lb  +  o-cao-cb  pab 


(11-17) 


from  ( II - 7),  where  pab  is  the  correlation  coefficient  between  the 

calculated  values  Lca  and  L>cb  for  Paths  a  and  b,  respectively.  Then, 

If?  Z 

P]_3  ”  o-^a  o"cb  t*3-  Woa  ^a  ^b  +  ^ta  wca  ^a  ^b 

+  ata  wCa  °ca  "-cb  Pab  "  “tl  woa  La  Lb 

-"tawcaLaLb}  =^pab  (11-18) 


p^is,  of  course,  identically  equal  to  Pab« 
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Substituting  now  (11-13), 
°tb  =  +  ac l  +  ^cb  ‘ 


(H-18), 

2<rta 


and  pal5  into  ( II- 6) , 
^ta 


ta 


+  2crta  acb -  Pab  "  2  aca  ^cb  Pab 


ca 


^  2  ,  „  2 
"ca  +  ^cb 


17 ta  +  ^caVbPab  '  l)  (H-19) 

ca 


If  the  two  paths,  Path  a  and  Path  b,  are  so  different  that  the 
correlation  coefficient  between  the  calculated  values  may  be  set  equal 
to  zero,  (11-19)  reduces  to 


"tb  =  °”ca  +  7cb  '  7ta  (II-20> 

In  this  case  no  improvement  for  Path  b  is  obtained  on  the  basis 

of  measurements  for  Path  a,  as  c rt£  will  always  be  larger  than  <r  The 

latter  is  the  prediction  uncertainty  of  the  calculated  value  for  Path  b 

without  any  measurements,  and  will  in  this  case  be  more  advantageous 

to  use  than  cr^.  The  improved  prediction  uncertainty  crta  for  Path  a 

will,  of  course,  be  always  smaller  than  o"ca,  based  only  on  calculated 

values,  and  the  right-hand  side  of  (11-20)  will  always  be  larger  than 
2 

^cb* 

If  the  Paths  a  and  b  are  completely  identical,  p  ^  =  1,  and 
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Then  (11-19)  reduces  to  the  identity 


(11-21) 


The  case  of  most  interest  in  practical  applications  is  the  one 

where  the  paths  are  nearly  equal  (i.e.  ,  in  the  same  general  area  of 

similar  meteorological  or  terrain  characteristics,  with  perhaps  one 

of  the  terminals  identical).  In  this  case  the  prediction  uncertainties 

cr  ^  and  cr  ^  are  equal^/ 
ca  cb  n 


(11-22) 


and  (11-19)  is  transformed  to 


"t l  =  *ta<2Pab  -  ‘>  +2lrc(1  -P>  <n-23> 


Dividing  by  cr^,  and  denoting  the  prediction  improvement  ratio 
for  Path  b  cr^/o-^  as  r^,  r^  may  be  obtained  as  a  function  of  the  pre¬ 
diction  improvement  ratio  for  Path  a  (for  which  the  measurements  are 
available)  cr^/tr^  ~  and  the  correlation  coefficient  Pa^' 


r2  =k2(2pab  -  1)+  2(1  -  pab)  (11-24) 


Curves  for  r^  have  been  calculated  and  are  shown  on  Fig.  II  - 1 . 

2/ A.  P.  Barsis,  K.  A.  Norton,  and  P.  L.  Rice,  "Predicting  the  per¬ 
formance  of  long  distance  communication  circuits,  "  NBS  Report  6032, 
December  1958 
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Only  values  of  r  <1  represent  an  improvement  in  the  prediction  for 
Path  b  <  tr£).  Fig.  II -1  shows  that  pa^  has  to  be  at  least  0.  5,  and 
the  improvement  for  Path  b,  of  course,  depends  on  the  degree  of  im¬ 
provement  =  o't^/o'^  obtained  by  measurements  for  Path  a.  For 
k  =1,  drops  out  of  (11-24)  and  r  =1  for  any  value  of  pa^. 


Figure  II  -  I 
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